Deformed two center shell model 
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A highly speciahzed two-center shell model has been developed accounting for the splitting of 
a deformed parent nucleus into two ellipsoidaly deformed fragments. The potential is based on 
deformed oscillator wells in direct correspondance with the shape change of the nuclear system. For 
the first time a potential responsible for the necking part between the fragments is introduced on 
potential theory basis. As a direct consequence, spin-orbit Is and operators are calculated as shape 
dependent. Level scheme evolution along the fission path for pairs of ellipsoidaly deformed fragments 
is calculated. The Strutinsky method yields the shell corrections for different mass asymmetries from 
the superheavy nucleus ^°^122 and ^^^Cf all along the splitting process. 



INTRODUCTION 



From more than thirty years on, two-center shell models are precious tools in the research of fusion, fission and 
cluster decay processes. The capability of producing transition level schemes during the splitting process from an initial 
parent nucleus towards two different fragments offers the opportunity of studying what are the possible microscopic 
changes through a certain fission or decay channel. 
. In 1961 a symmetric double center oscillator potential has been analyzed by Merzbacher 0. A strong link was 
emphasized with regard to diatomic molecules (like ammonia molecule) where the motion in the neighbourhood of 
^ I an equilibrium state is close to the one generated by a harmonic potential. This problem has been solved for two 
' identical spherical potential wells. 
\^ In the 1970's, J. Maruhn, W. Greiner and the Frankfurt school developed the asymmetric two-center shell model 

0- It was an important step for the study of mass asymmetry in binary nuclear fission. The model considered two 
, spherical asymmetric fragments subjected to two-oscillator type potential. A remarcable feature of this model was 
' its ability to describe two overlapping spherical nuclei fission shapes up to their separated configuration. Thus it was 
possible for microscopic effects to be held responsible as being specifically due to the two emerging and gradually 
separating level schemes. Basically it is under this form that two center shell models are widely used today in nuclear 
fusion, fission and cluster decay phenomena. 

Two center states in light nuclei are taken into account in connection with molecular orbitals in a two-center shell 



^ model picture 0]. Also, application of the two-center shell model on the study of ^^0-|-^^C reaction has been published 



with emphasis on nuclear Landau-Zenner effect. 

Nucleus-nucleus collisions have also been considered in the framework of the two-center shell model to obtain the 
' two-center level diagrams and single particle corrections for asymmetric systems such as ^^'^'^O _|_24,25 j^j^g gj 
^ . The importance of an adequate description of cold fission, cluster radioactivities and alpha decay in terms of an 
' asymmetric and deformed single particle shell model with more realistic shapes during fission and fusion processes 
was repeteadly stressed 

Recently, attempts have been made to use the two-center shell model in synthesis and decay of superheavy nuclei. It 
is a suitable theoretical model to study the microscopic effects on possible projectile-target combinations in their way 
from two different quantum systems to one. As the synthesized nucleus is heavier, increasing Coulombian repulsion 
lowers the macroscopic potential barrier almost to complete dissapearance. The only way a superheavy element can 
survive is due to the shell effects |^]. A two-center shell model is essential for the description of fusion and fission 
of superheavy elements. It shows that the shell structure of the two participating target-projectile nuclei is visible 
far beyond the barrier into the fusioning nucleus and is crucial in choosing the most favorable pair approaching each 
other through cold fusion valleys. Transition behaviour of the two partner shells will provide shell corrections which 
lower the fusion barrier to be overcome, as compared to neighbouring projectile-target combinations. 

Up to now, all the variants of two-center shell models used spherical nuclei. I shall consider now the motivation of 
this work. In any process which implies a pass from two quantum systems to one (fusion) or the other way around 
(fission) there are certain situations when one or both fission fragments or fusion partners are deformed. Such a 
reaction could yield cold energy valleys due to deformed shell structure of the participants. Fragment deformations 
are properly accounted in this work. A deformed two-center shell model (DTCSM) is proposed, where the main part 
of the potential consists of two ellipsoidaly deformed Nilsson type oscillators for axially symmetric shapes. Any change 
in the nuclear surface shape is reflected in a corresponding modification of the four oscillator frequencies along the 
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symmetry axis and perpendicular to it. 

It is also well established that, especially in fission and cluster decay, a necking region builds its way between the 
fragments, smoothly linking the two ellipsoids one to the other. This work associates for the first time a microscopic 
potential to a spherically matching neck region of the nuclear shape. The neck potential is constructed in such a 
way that it takes the same value on the necking region surface as on the ellipsoidal region surfaces of the fragments. 
Equipotentiality is thus respected on the nuclear surface. 

The usual spin-orbit and squared angular momentum operators are calculated with help of the potential-dependent 
formulae Is — (VF x p)s, and P = (V^ x p)^ The potential in this model follows exactly the nuclear shape, 
and so do the Is and 1^ operators. Finally, the levels diagram with respect to the elongation, neck parameter and 
fragment deformations are used to calculate the shell correction by means of the Strutinsky method. 

Calculations are presented as an example for diff'erent mass asymmetries, for the superheavy nucleus ■^"^122 and for 
252q£ rpj^g dependence on various neck parameter values and on different pairs of fragments is discussed. 



SHAPES 



Fig. 1 shows the main geometrical parameters defining the axially symmetric shape family DTCSM is dealing 
with. Two ellipsoids (the deformed fragments) with semiaxes ai, bi and 02, 62 are, at a certain moment, separated 
at a distance R between the two centers Oi and O2. A sphere centered in O3 with radius R3 is rolling around the 
symmetry axis, being tangent all the time to the two ellipsoids. The necking region, between the two tangent points, 
is generated in this way. Thus we have five independent parameters to design the deformation space: two fragment 
shape asymmetries xi = bi/ai, X2 = ^2/02 (if and 02 are given as the correspondent of [^2 for every Ai and A2, 
the other semiaxes are calculated from the total volume conservation condition), mass asymmetry A1/A2, the neck 
radius R3 and the distance between centers R. Obviously, this set is available for every parent nucleus A, Z with its 
own X = b/a. 

A few shape sequences obtained by varying two of the parameters, R3 and R for the same parent and the same 
mass asymmetry A1/A2 are depicted in Fig. 2. Every row starts with the same ellipsoidal parent on the left hand 
side of the figure. As the distance between centers increases, the two deformed fragments with fixed xi and X2 
shape asymmetries begin to separate one from the other. Variation with the neck radius R3 is noticeable on vertical 
direction. Shape sequences with very small necking region are shown in the upper row (for zero neck radius, R3 = 
fm, we get compact shapes suitable for fusion reactions), passing through intermediary neck radii, comparable to the 
magnitude of the ellipsoids semiaxes (second and third row), down to the last row where large neck radius generates 
very elongated shapes. 

This kind of configurations will be microscopically treated along the variation of the deformation space parameters. 



THE POTENTIAL 



The equations for shape surfaces described in the previous section can be written in cylindrical coordinates (due to 
axyal symmetry) as: 



p{z) = <^ pg{z) = P3 - [Rl - (Z - Z3)2]l/2^ Z,i < Z < Z,2 (1) 
yp2{z)= [hl~X2{z-RYY'^, Zc2<Z<R+a2 

where the origin is placed in the center of the heavy fragment Oi. Neck sphere center coordinates are (23,^3), and 
Zci and Zc2 are the two tangent points of the neck sphere with the two ellipsoids. 

The oscillator potential correspondig to these two-center shapes must have the same value on the nuclear surface. 
For spheres , for example we have: 



^0 - ^ (2j 

where Ri is the radii of a nucleus with atomic mass Ai. Since huji = ilA'^/"^ and Ri — tqA^^^ (where tq ~ 1.16) 
then Vo ~ 54.5 MeV. If we write these simple relations for the surface of ellipsoidal shapes: 
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-mo<6? = Vo (3) 



the frequencies ujz. , ujp. are defined along the symmetry axis and respectively perpendicular to it, as functions of the 
two ellipse semiaxes. 

For an arbitrary origin, placed on the symmetry axis, the ellipsoids surface equations read: 

p2 ^ jz + Zl)^ ^ ^ 



2Vo ' 2Vo 



„2 ^\2 



2Va ' 2Vr 



where zi and Z2 are the absolute values of each of the two centers coordinates. Now the two oscillator potential 
expressions for deformed fragments come straightforward; 



1 2 2 1 

1 o o 1 



yiiP,z) = -moujl^p^ + -moujl^^z + zif 



y2{p,z) = ^moojiy + -moojl^iz - Z2f (5) 

What we have left to establish is the necking region potential, Vg{p,z). The force that keeps nucleons confined 

within the ellipsoid is an attractive type one, radially inward {Fpiupsoid ^ The same reason leads us to the 

hypothesis that, if nucleons are confined within the concave necking region, which is geometrically inversed to the 
ellipsoid convex surface with respect to the centers of the fragments, a rejective force is needed, radially outward: 

F(r) = -/3f (6) 

hence a force oriented from outside the nuclear shape toward the surface. Then the corresponding potential is related 
to the expression: 

/ Fdf = <fi{r) - (p{0) (7) 
Jo 

or 

/ Fdf= [ -l3rdf=-^ (8) 
Jo Jo 2 

Consequently, ip{r) = — is a potential form generating the rejective force. In this way, the rejective neck potential, 
defined up to a constant, must look like: 

Vg^{r)=V, + ^{r) (9) 
where 14 is a constant to be determined. On the nuclear surface Sg, we have again: 

Vgi\s=Vc + ip{r)\s=Vo (10) 

or 
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= ^0 " Vir) Is, (11) 
For the sake of consistency, the rejective force is also considered of osciUator type: 

F{r) = -moLulf (12) 

therefore 

V>{r) = ^ (13) 

where the frequency LUg has to be found. Since the potential must follow the geometrical shape, at the neck region 
the function ip{r) — (p{p, z) reads: 

^{p, z) = -^[(p - P3)' + z,?] (14) 
and is centered in the middle of the neck sphere 03(^3, Z3). On the neck surface then, where (p,z) £ ^g, we have: 

^{P, z) -^Rl - Vgi l5,= ^0 (15) 
Then Vc — Vq — Lp(r) \sg— 2Vb and the total neck potential from outside the shape down to the surface is: 

Vg^{r) = 2Vo - ^[{p - p,f + {z - z,f] (16) 

and the neck frequency is directly related to the neck radius by Eq. |l^. Vgi reach its maximum at the center of the 
neck sphere {p — p3, z — 23), where Vgi = 2Vo, then is decreasing down to the surface value Vgi = Vq at a distance 
equal to the neck radius (p — ^3)^ + [z ^ z^Y — ^3 from O3. 

To complete the neck-dependent potential, there still remains the region inside the nuclear shape between the 
necking surface and the interior contours of the two ellipsoids (black colored, denoted by Vg2{p, z) in the upper part 
of Fig. 3). It can be observed that on the ellipsoids surface inside the shape the deformed oscillator potential has the 
same value as on the fragments surface, namely Vq. But on the surface of the shape within the necking region, the 
value is also Vq. Then, one concludes that inside the region volume between the neck surface and the two ellipsoids 
surfaces - pi{z) < p < Pg{z) and P2{z) < p < Pg{z) - the neck potential is constant. 



Vg2{p,z)^cst = Vo (17) 
Finally, the deformed oscillator potential part for the DTCSM reads: 

{Vi{p,z)^ imoWpjP^ + imow2^(z + zi)2 , 

Vgip^z)^ ( 2K,-[i™o-^(p-p3)^ + imoa;,^(z-Z3)^] ,.,1 ^^3^ 
t Vg2{p,Z) - Va ,Vg2 

V2{p,z)^ ^niQUjI^p^ + \moUjl^{z - Z2Y , U2 

where wi, Vgi^ Vg2 and V2 are the spatial regions where the corresponding potentials are acting. These regions are 
about to be revealed further on. 



MATCHING POTENTIAL SURFACES 



One of the most important issues for the two-center shell models is that the separation plane between the two 
fragments is not the one where the potentials equal eachother. As an example, for compact shapes (no neck, R3 = 0) 
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the points ( {ps,Zs} ) on the separation line (where the two elhpsac intersect) result in different values for the 
potentials, Vi{pg, Zs) ^ V2{ps, Zg). The sharp cusp between the two values makes calculation totally wrong. Attempts 
to solve the problem by a geometrical transition function, making the frequency ojp^{z) join LOp^{z) introduce an 
unpredictible approximation which is bigger as the mass asymmetry is larger; moreover, it does not reproduce the 
neck spherical shape and still, only the p— part of the potential is solved. 

In conclusion, one has to divide the Hilbert space such that the whole potential Vdtcsm{p, z) could smoothly and 
continuosly pass from one fragment to the neck and finally to the other fragment. 

First, let us take the deformed fragment oscillators into consideration and find out what are the spatial limits 
of continuity between the two potentials. The equation for such a region must be complacent with the matching 
condition: 

vMz) = y2{p,z) (19) 

where 



V2 = ^moLol^p^ + \mQUjl^{z - Z2f , 5'2(p, z) 

or 



(20) 



py + < + ZiY = + <(2 - Z2Y (21) 



If one translates the 2;— coordinate: 



z = z' -k (22) 
and separate the terms with respect to the powers of p and z' , we get: 

« - O + K - Oz" - 2M< - ^ly + K - 

+2{ivlzi + L0%Z2)z' - 2{ulzi + Lol^Z2)k + ivlzf - u^zl = (23) 

We notice that: 

zi+Z2 = R (24) 

and after some simple calculations one obtains: 

+ ^f^ = l (25) 



, ,2 , ,2 p2 ^ , ,2 , ,2 p2 



which defines a so called unique center surface, here namely an ellipsoid with semiaxes: 



^ uj^^ 

Since the origin of this surface is at (p = 0, z' = 0), its center with respect to the arbitrary origin (against which 
the potentials are expressed) is on the symmetry axis at: 
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a\z2 



(27) 



What wc have obtained are the coordinates (zomjO) and the semiaxes (a„j.6„j) of a matching potential eUipsoid 
(MPE). On its surface Vi{p,z) = V2{p,z). The heavy deformed fragment Oi is acting through Vi{p,z) outside the 
MPE. The hght fragment part which emerges from the parent is contained inside MPE (the shaded region in Fig. 3). 
Hence its action on the Hilbert space, V2(p, z), goes all over the matching ellipsoid volume. 

The same demonstration is valid for Vi{p, z) = V^i(p, z), and V2{p, z) = Vg2{p, z), as far as the neck potential shares 
its action space with the fragments. First condition: 



or 



yields: 



Vi{p, z) = Vgi{p,z) 



^moJl^p^ + \'m.oujl^{z + zif = 2Vo - [^moWg(p - ps)^ + ^moU!g{z - zsf 



(28) 
(29) 



(moWpj + moUjg)p'^ + {moLul^ + moWg)z^ - 2moLUgP3p + 
2{mou'^^zi - mo 23)2: + moLol^zl + moajg{pl + zj) - Wq = 

This equation is of the type: 



(30) 



aiip'^ + a22z'^ + 2aip + 2a2Z + a — (31) 
To describe a unique center surface, the normalized determinant of the equation: 

A a? Oo 

-=a-^-^ 32 
S ail (122 

must be positive, which is the casefor our situation. Then the center of the surface has the coordinates: 



Pomi — 



After some simple calculations one obtains: 



POml 



ai2 ai 

a22 0,2 



This surface is again an ellipsoid with semiaxes: 



an oi 
021 02 







ZQvaX — 






an 


ai2 




an 


«12 


a2i 


a22 






ail 



(33) 



ZOml 



U)gZ3 - i^l^Zl 



(34) 



0-m\ 



A 1 



1/2 



O22 



A 1 

8 ail 



1/2 



(35) 



and represents the matching potential ellipsoid (MPEl) between Vi{p^ z) and V^i(p, z). On its surface, the two values 
equal each other. Obviously, the tangent point between this first matching potential ellipsoid and the heavy fragment 
shape is the same as the tangent point between the neck sphere and the heavy fragment. In the interior of MPEl, 
the DTCSM potential is V^i(/0, z). Following the same arguments for: 
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V2{p,z) = V,i[p,z) (36) 

one finds: 

P0m2 — ^ ; 5- 20m2 — 5 j ^ (,0 ' j 

for the center of the second neck matching potential elhpsoid (MPE2), and the same expressions with V2(/3, z) corre- 
sponding values instead of Vi hold for the semiaxes 0^2 and 6^2 ■ 

In the upper part of Fig. 3, the three matching ellipsoids - the large shaded one for Vi = V2 (MPE) and the two 
grey ones for Vi — Vgi (MPEl) and V2 = Vgi (MPE2) respectively - are drawn. In this way DTCSM potential spans 
the whole Hilbert space as follows: within MPE it is V2(p, z) which is active; within the grey area formed by the two 
neck potential matching ellipsoids, MPEl and MPE2, we have Vgi(p, z); in the interior of the nuclear shape, between 
the neck surface and the fragment ellipsoids surfaces (tiny black region) the potential is Vg2{p,z). The rest of the 
space is under Vi{p, z) action. 

What is achieved in this section is spanning the Hilbert space by DTCSM potential, passing continuously from the 
potential generated by the heavy fragment Vi to the lighter one V2 and allowing also for the neck to be smoothly 
shaped by Vg. There are no sharp cusps in the potential values. Since we work in axial symmetry, the two neck 
matching ellipsoids MPEl and MPE2 are in fact rolling around the symmetry axis. Each of them shapes a thorns 
around the necking region, where Vgi is the DTCSM potential. 

On the lower part of Fig. 3 geometrical points which will be used for further calculations are marqued: Pmis{z) 
and Pmij{z) are MPEl surface functions, Pm2s{z) and Pm2j{z) hold for MPE2, and Pm{z) is the surface of MPE. 
Other points of interest are: Zd and Zc2 - the two tangent points of the neck and MPEl and MPE2 with the fragment 
ellipsoids, - the intersection of MPE with the symmetry axis, Zxu and Zx2r - the left margin of MPEl and the right 
margin of MPE2, Zg - the crossing point between the two fragment ellipsoids, Zmint - upper crossing point between 
the three matching ellipsoids 

Fig. 4 shows how the three potentials work together. The values of Vi(p, z), V2{Ptz) and Vgi{p^z) are taken at 
p = pi(zci) (upper part of the figure), and p = P2{zc2), hence at the tangent points of the neck sphere. Evolution 
is shown as the distance R (or z-variable) increases. The values are taken at a matching surface point, thus the 
potentials are tangent all the time. One can observe how Vi and V2 wells are separating with increasing R. The 
matching surfaces are composed of the whole set of tangent points at every p distance from the symmetry axis. 

The way two of the potentials, Vi (p, z) and Vgi (p, z), behave along the O1O3 line between their two centers is given 
in Fig. 5. The crossing points of the curves correspond to Zd and its opposite along MPEl. Inside the neck sphere, 
hence inside MPEl too, Vgi > Vi. 

DIAGONALIZATION BASIS 

This section will browse succintly the main steps to be taken in order to find a normalized functions set, since 
details about this issue have already been published § . 

The new feature in the DTCSM eigenvalue problem is the shape-dependency of the potentials: Vdtcsm, the spin- 
orbit term Vos and the V02 term ($7 is preferred here as a notation since the angular momentum operator, as being 
shape dependent, will be different from the usual 1). The total Hamiltonian: 

^2 

H DTCSM = -7; — A + Vdtcsm{p, z) + Vhs + Vb^ (38) 
ZniQ 

is obviously not separable. A basis is needed and diagonalization of oscillator potential differences and of angular 
momentum dependent operators has to be performed. 

A separable Hamiltonian is obtained if one takes — = cji, with no Is and F terms, hence a potential like: 

1/W(p z) = / ^i"*^*^^'^) = ^moUjfp^ + ^moUjf{z + Zi)'^ ,z<0 ^^g^ 
)vj'^\p,z)= ^uiqujIp^ + \m.QUjl{z - Z2Y ,z>0 
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This is an appropiate two-center potential for a sphere {z < ) intersected with a vertical spheroid. The origin 
(z ~ 0) is the intersection plane. As a result of variable separation, three known differential equations are obtained 
for harmonic functions, Laguerre polynomial and Hermite function dependent solutions 

Using continuity conditions for z-dependent functions and their derrivatives at z = 0, the normalization condition 
and equalizing the energy on the symmetry axis Ez = Tiwiii'i + 0.5) = 1iuj2{v2 + 0.5), z-quantum numbers and 
normalization constants for the z-dependent eigenfunctions are calculated. The ^ and p-dependent functions are 
straitforward calculated from the differential equations. The final result is: 

^m('/') = ^= exp (im^) 



27r 



Z,{z) = 



a,cxp -fiil£+£lL H,,[-ai{z + zi)] , z<0 



C^2 exp 



Hi,.^[a2{z - Z2)] , z > 



We notice that: where T is the gamma function, L™is the m-order Laguerre polynomial, Ci and C2 the normalization 
constants, vi, V2 the quantum numbers along the symmetry axis, and Hi, is the Hermite function. We now have 
a basis for further calculations. The eigenvalues for the diagonalized Hamiltonian with the potential V'^'^^ are the 
oscillator energy levels for sphere+sphcroid system: 

Efl = hujx i2np+ I m I +1) + huj^, {1^1 + 0.5) (41) 

At this point we have a useful basis for the calculation of Hdtcsm matrix elements. 

Hdtcsm operators 

In order to obtain the DTCSM energy levels, the matrix of the non-diagonal elements should be constructed. Then, 
adding E^ffJc as diagonal terms, after diagonalization we obtain the eigenvalues. 

Vdtcsm{p,z) - terms 

This subsection is devoted to build the operators for the necked-in deformed two-center oscillator levels calculation. 
We start wiht the operators leading to the appropiate energies for Viip, z). The difference Vi{p, z) - V'''^\p, z) on the 
volume Vi and added to E'^'^\np, i>,m) produce the suitable matrix elements for Vi: 

AVi = AVi{ ) + AVi{ 0) -I- Al/i(0 ) + AVi{zQra + ajn,oo) (42) 

-AVi{Zj:ii Zci;pmlj ^ Prnls) - AVi{Zcl ^ Zs\ pi Pmls) " AVi{Zs ^ ZminU Pm Pvals) 

The positive terms account for the influence of Vi(p, z) as if Vi would act alone everywhere but MPE. The negative 
terms are under Vg{p,z) action. In parenthesis are the geometrical limits for every term. Note that for those which 
have also p-dependent limits, p — p{z). The limits can be easily tracked in Fig. 3. What is left after the subtraction 
is the difference from Vi{p, z) to V'^'^\p, z) within the vi volume. To ease calculations, one uses as a notation: 

0) 4- AFi(0 ) -f AVi(zo,„ -f a,„,oo) (43) 

These pozitive potential differences are now written in terms of the Heavyside step function; this will result in 
determining the exact limits where Vi{p, z) is acting and, consequently, where it decides the oscillator energy values : 

^^1^^ = i(moWp^ - mocj^){[l - 6l(z- z,„0] + 6i[z- (zo™ + a™)]}p2 
+ i{(mowf^ - moWi)(z + zi)^[l - 9{z - Zim)] 

+ [moujl_^{z + zif - triQUjliz - Z2)^]6'[z - (zo,„ + a™)]}^„^,„p (44) 
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+ -(moWp^ - mnu;l){[6{z - Zi,n) - 9{z)] 

+ [e{z) - 9[z - (zo„, + a^)]]}p^9{p - p„.(z)) 

+ -{(mow^^ - moWi)(z + zif[0{z - z^m) - 6{z)] 

+ [mowfj(z + zif - moujliz - Z2Y][0{z) - 9[z - (zom + ayn)]]}6{p - Pra{z)) 
For the negative differences we have: 



/\Vi{z^ ^ Zf] p,{z) ^ Pf{z)) = 

imo< [6{z ~ z,)6{z ~ - p,(z)) - 6[p - pf{z))] 

+\m^Lol{z + z,f[e{z - z,) - e{z - zf)][e{p - p,{z)) - e{p - pf{z))] 



(45) 



where Zi, Zf, pi{z) and Pf{z) tafce the corresponding values and hmit functions from Eq. 45. 

To obtain the terms where the hght emerged fragment potential V2{p,z) acts, one must subtract V^'^^{p,z) from 
V2{p, z) within the V2 volume. The result is: 



Let us denote: 



AV2 = 

AV2izim 0; Pm{z)) + AV2iO zo,n + a„ ; p„i(z)) 

~AV2{Zs Zc2\Pg{z) ~ Pm{z)) ~ AV2iZmint ^ Z^2r\ Pm2]{z) P,n2s{z)) (46) 



AV^'P^ = AV2{z^„, - 0; - p^{z)) + Al/2(0 ^ zom + a™; - p™(z)) (47) 

Then we have: 



-{moujp^ ~ mQUjl){0[z - (zom + a™)]}p^[l - 9{p - Pm(z))] 
+ i{[mowf^(z - Z2)^ - ruQUjKz + zi)^][6l(z - z™) - e{z)] 

+ (mow^^ - mow^)(z - Z2)'^{e{z) - 9[z - (zq™ + a,„)]}}[l - 9{p - Pra{z))] (48) 

for the positive part; for the negative necking region which has to be subtracted from V2, the expression for any 
AV2(zi ^ Zf;pi{z) ~ Pf{z)) is: 



AV2iz,~Zf;p,^Pfiz))= (49) 
1 

—1 

+ imo<(z - Z2n9{z - Zi) - 9(z - Zf)][9{p - pdz)) - 9{p - p;(z))] (50) 



-mo<[0(z - Zi) - 9{z - zf)]p^[9{p - p,{z) - 9{p - pf{z))] 



using the same convention for Zj, zj, pi and Pf{z) as for AT^i. 

At this point, the neck potential Vg{p, z) is filling its volume Vg without any subtraction: 



AVg = Vg 

With the same notation as in Eq. we have for the neck operators: 



(51) 



10 



Vg = Vgl{Vgl)+Vg2{Vg2) 



(52) 



Now every of them will be taken separately: 



Vgl — Vgi{ZxU -T- Zci;pmlj -j- Pmls) + Vgi{Zcl -T- Zmint'i Pg Pmls) + Vgi{Zm,int Zxir'i Pm2j Pm2s) (53) 



where again: 



and 



Vg2 = VgliZcX ^ Za\p\{z) ^ pg{z)) + Vg2{,Zs ^ Zc2\ P^iz) ^ Pg(z)) 



Vgi{zi^ Zf;pi^ Pf) = 

2V^[6{z - z,) - 6{z - Zf)][e(p - pdz)) - 6{p ~ pf{z))] 

-^-moujlnz - Zi) - e{z - zi)\{p - p^f[e{p - pi{z)) - e{p - pf{z))\ 

-\mou^l{z - z^f[e{z - Zi) - e{z - zf)][e{p - pi{z)) - e{p - pf{z))] 



Vg2{zi - zf, Pi - pf) = Vamz - Zi) - 6{z - Zf)][e{p - pi{z)) - 6{p - pf{z))]} 



(54) 



(55) 



(56) 



The final analytical expressions for the matrix elements corresponding to the oscillator potential of DTCSM are given 
in Appendix A. 

Spin-orbit Is and 1^ operators 

Special care will be devoted to the treatment of Vis and V\2 matrix elements. Because of the dependence of the 
angular momentum term on different space regions ( ~ W x p), hence on different mass regions when asymmetry 
A1/A2 comes in, the anticommutator is used to assure hermicity for the operators: 



Vis 



mowoi 

n 

TO0W02 



-Ki{p,z),{VV X p)s 

K2{P,Z),{VV Xp)s] 



, Ai — region 
, A2 — region 



(57) 



and 



V12 = { 



-^Kipi{p,z),{VV X p) 

-2-^K2M2(p, ^),(VFxp) = 

. '"'0^02 



, Ai — region 
, A2 — region 



(58) 



Basically the same treatment as for oscillator terms is valid in this case. The Ki{p, z) and /ii(p, z) are the strength 
function parameters for the V\{p, z) region, whereas K2{p, z) and /U2(p, z) are active for the V2{p, z) one. 

The new feature here is the use of the anticommutator for the operators containing Hcavyside function combinations, 
to confine the action to wi, V2 and Vg] these function combinations are exactly the ones which have been used for A^i, 
AV2 and AV^ terms. 

Since the creation and anihilation angular momentum operators become shape-frequency dependent, notations will 
be changed as follows: 
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so that 



1 



Then for the potentials of the spin-orbit term to be diagonahzed one reads: 



(59) 



where: 



Vhs = Vhs(wi) + VnsM + Viisivg) 



(60) 



Vhs(fl) 



mot^oi 



-Ki{ns, (vi)} 



(61) 



mouJo2 



-K2{rjs, (W2)} 



(62) 



V^s{vg) = —Ki{ns, (z;«)} —K2{ns, («<,)(2)} (63) 

niQUJoi " moujo2 

where Wg^^ and Wg^'' are the neck matching eUipsoids volumes on Ai and A2 side, respectively. 

On the other hand, the operators depend on the region they exert themselves through the potentials, according to 
Eq. pn. Thus we have: 



n+{vi) = -e'f 



dVi{p,z) d dVi{p,z) d idVi[p,z) d 
dp dz dz dp p dz dkp 



d 



d 



a 



'^o^piP'gZ ~ "^o^^ziiz + zi)— - -moUj^_^{z + zi) — 



dp p 



dip 



(64) 



dVi{p,z) d _ dVijp.z) d ^ i dVx{p,z) d 
dz dp p dz dip 



dp dz 
d 



d 



d 



"^ot^Pi^^ - mow.^iz + zi)— + -moiJ,^{z + zi) — 



dp p 



dip 



(65) 



with 



^z{vi) = 



i dVi d 
p dp dip 
d 



dip 



(66) 



rjs(ui) = -(r2+(ui)s +fl s+) + flzivi)sz 



(67) 



The same aplies for (U2) and (vg) spin-orbit operators terms, and the expressions are: 
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mou) [z - Z2) 



dp 
d 



i 2 I ^ d 
-mou,,^{z-Z2)g^ 



2 d 
d 

^0^P2p-g^ - '^0<(^ - Z2)-^ + ^mo<(^ - Z2)-g^ 



(68) 



and: 



moulip - pz)— - mowl{z - zo.)— - -moul{z - z^) — 



2 P-Pi 9 

Q O 



(69) 



and fis(ug2) = 0, since Vg{vg2) = Vq = est. With the help of these operators, what is to be calculated now for the 
spin-orbit shape-dependent matrix elements reads, for [vi): 



n 



mowoi 



Ki{^s{v{), (wi)} 



(70) 



where: 



{ns(wi), [1 - e{z - Zim)]} + {ns(wi), e[z - {Zom + Om)]} 
+ {ris(wi), [e{z - Zim) - 0[z - (20m + am)]]0{p - Prn{z))} 
-{fts{vi), [e{z - Z^ii) - e{z - Zcl)][9{p - Pmlj{z)) - 0{p - pralsiz))]} 

-{ns{vi), [9{z - zei) - 0{z - z,)][e{p - pi(z)) - 9{p - PmUz))]} 

-{ns(?;i), [61(2; - Zs) - 0{Z - ZmintWip - Pm{z)) - 0{p - Pmu{z)]} 



(71) 



For the V2(p, z) controlled region we have: 



Vns{v2) = — 



n 



mou)o2 



K2{^S{V2), {V2)} 



(72) 



where: 



{f2s(«2),(«2)} = 

{r2s(w2), [e{z - - e[z - (zom + a„0]][l - 0{p- p„i{z))]} 

-{ns{v2), [e{z - z,) - e{z - z,2)mp - p2{z)) - e{p - pg{z))]} 

-{ns{V2), [6{Z - Zgm) - 6{Z - ZmintWip " Pg{z)) " ^{P " Pm{z))]} 

-{ns(W2), [0{Z - Zmint) - 0{z - ^x2r)][^(P " Pm2j{z)) - 6{p - pm2s{z))]} (73) 

Finally, the neck potential dependent spin-orbit interaction looks like: 

Vns{vg) = —K,{ils{vg,), (««)} - -J^,2{n^{vg,), {vf;;)} (74) 

mowoi moWo2 

where: 
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{ns{v,,),{vl\^)} = 

{ns{vgi), [0{z - z^ii) - e{z - Zci)][9{p - p„aj{z)) - 9{p - Pmis{z))]} 

+ {ns{Vgi), [e{z - Zcl) - 9(z - ZmintW{p~ Pg{z)) - e{p - Pmls{z))]} (75) 



{JlsKi),(«fi^)} = 

{r2s(Ugi), [e{z -~ Zraint) " 0(z - Z^2r)][d{p - Pm2]{z)) - 0{p - Pm2s{z))\} (76) 

For the P — > fJ^ term the usual expression works: 

qF- + ct-n+) + rLl (77) 

and the operator for the whole space is: 

Vn-=-\\^^^^^P^{p,z),n'\ ,i = l,2 (78) 

Strength coefRcients z) and K2P2{p, z) are the same {p, z) - dependent functions as Ki(p, z) and K2{p, z) in 

the spin-orbit term. Making use of the unity operator / = ^„ |")("|, many expressions can be derived to obtain an 
easy workable decomposition of the lengthy matrix elements of Vfi2 . The chosen one follows, in order to make use of 
the already calculated expression for Jls term: 

('\{fip,z),n+^i-}\) = 

^[2('|{/(p,z),f2+}|") + ('|(f2+/(p,z))|")]("|f^-|> + ('|f^+|">("l(J^-/(P,^))l) (79) 



{'\{fip,z),n-n+}\) = 

^[2('|{/(p,z),Jl-}r) + ('|(fl-/(p,z))r)]("|f^+|) + ('|f^-|")("l(f^+/(p,^))l) (80) 

// 

The needed expressions for the fls and $7^ matrix elements are detailed in Appendix B. 
Now the total matrix elements for DTCSM can be calculated as: 

{t\DTCSM\j) = Eiiliup, |m|,/i) + {i\AVi\j) + {i\AV2\j) + {t\Vg\j) + {i\Vns\j) + m^-lj) (81) 

LEVEL SCHEMES AND SHELL EFFECTS 

A ten shell scheme matrix elements have been calculated, which means a 440 x 440 matrix. After proper diago- 
nalization, the levels are obtained. 

First, DTCSM spectra are computed for the superheavy fission reaction channel ■^°^122 ^^^^ W -{-^'^^ Cd, with 
the nuclei deformations X122 — 0.9, xw — 1- f^nd xcd — 0.83. Semiaxes and bi are calculated from corresponding 
deformation parameter (32 for every fragment, and using the volume conservation condition {(32 ^ =0., /?2 =0.135 
). The reduced distance between centers is chosen to represent the elongation parameter of the shape : (R — 
Ri)/ {Rf — Ri), where Ri is the distance between centers when the light emerging fragment is completely embeded in 
the parent nucleus, and Rf — ai + a2 + 2R3 represents the final distance between centers, when the neck sphere is 
aligned with the fragments. 

A direct consequence of the variation of the neck parameter i?3 on microscopic behaviour of a fission process is 
depicted in Fig. 6. Here the shell corections are drawn for the five R3 - parameter values, as a function of the distance 
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between centers. They are calculated with the Strutinsky method. A shallow minimum at i? « 2.5 fm suggest a 
small deformation of the ^^^122 ground state. First bump around i? « 5 - 6 fm shows up in every case. Second 
bump is also there, but its position changes in i?-value and in height (between R —12 fm and R —20 fm). The rather 
deep minimum (« -5 McV) for large neck radii (apparent on the R^ = 10 fm curve) is probably not manifested in 
306122 fission, since at this distance between centers, the system is out of the fission barrier due to strong Coulombian 
repulsion. 

Second example takes into account the mass asymmetry parameter. Again is the fission of the superheavy ■^^^122 
through two channels this time: 

306122 ^ i98^+i°8Cd 

The two corresponding level schemes are on the upper side of Fig. 7 {xGd =0.72 for ^2'^ = 0.243 and xce — 0.71 for 
(32'^ =0.261). They have been chosen because of the mass symmetry of the second reaction against an asymmetric 
one. Calculations in this case preserved R3 =4 fm. Shells are more clearly visible at {R — Ri)/{Rf — Ri)=l in the 
symmetric case, since asymptotically the levels are practically coincident for the two quasi-equal fragments. The 
lower part of the figure represents the shell corrections calculated with the above level schemes as an input. The first 
bump is shifted towards smaller R for symmetry. A second minimum as deep as almost the ground state along the 
symmetric channel is an indication of a possible shape isomer. Second maximum appears only for the asymmetric 
reaction. Asymptotically there is a AE «4 MeV difference in between the two channels, suggesting a possible fission 
path along the mass asymmetry degree of freedom. 

The last studied case is the fission of ^^^C/ {xcf =0.73 for (32 =0.236) through two favorized channels in Fig. 8: 

These two reactions have been chosen because they lie within the maxima of the fission fragment mass distribution of 
252c/ [jll]. The fragment deformation parameters considered in calculations are xSm=0-68 for /9f™=0.290; XKr—0-73 
for /3f '■=0.228; XBa=0.8 for /3f °=0.164; xmo=0.65 for ^^^"=0.333. Level schemes are rather likely in the first part of 
the splitting, a fact that can be observed from the small diferences in the shell corrections up to R «7 fm. Beyond this 
distance between centers, the eff'ect of the wells separation becomes more and more important, fragment individuality 
imposing the difference. If only for the shell corrections, S'm+^'^ Kr seem to be a more probable fission channel, but 
it beyond the purpose of this work to do a fission process analysis. Separation occurs around 18 fm, for approximately 
both of the reactions. As it is already known, lower energy shells of the fragments are the first to unveil asymptotically, 
whereas the higher states of the two systems still interact. This is an expected result since the potentials separate 
starting with the bottom of the wells, while the upper part of the potentials still overlap. 



CONCLUSIONS 



The deformed two-center shell model presented in this work describes the evolution of single particle levels from 
one parent potential well to the two fragments ones. 

The introduction of fragment ellipsoidal deformation in the two-center oscillator wells and further on, in the spin- 
orbit and l2 operators enables a more realistic calculation of the two interacting quantum systems. This new form 
of spin-orbit ils and fl^ allows for the angular momentum dependent operators to follow the exact sequence of 
shapes throughout the splitting process. In such a way ellipsoidal degrees of freedom are considered within spin-orbit 
interaction. 

The new necking-in dependent microscopic potential results in considering the neck degree of freedom into the level 
scheme calculation. As has been shown, the last part of the fision process is influenced by the difference in neck 
radii. This fact has important consequences on the potential barriers and flssion paths on a potential energy surface 
calculated within more deformation degrees of freedom, necking included. 

A new way of treating two partially overlapping nuclei has been introduced by making use of matching potential 
surfaces. It is only in this way that continuity between different regions of potential influence is assured. The potential 
pass smoothly from Vi to Vg , to V2 and so on without any cusp in its value and without the introduction of arbitrary 
geometric transition functions. For zero neck radius one obtains the fusion-like type of shapes. 
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All these facts make the presented model suitable for the study of fission channels and cluster decay phenomena 
and calculation of potential energy surfaces in fission. 

Also the model provides the analysis of possible deformed target-projectile doors toward fusion processes in heavy 
ion reactions. In nucleus-nucleus collisions, calculations of two fusioning deformed partners yield the single particle 
spectra for the important (yet avoided up to now) region of overlapping shapes. 



DTCSM OSCILLATOR MATRIX ELEMENTS 



Formulas for the two deformed oscillators with necking- in are presented to enable ('|AVi|), ('IAV2I) and ('|AVg|) 
matrix elements to be calculated. The spin-dependent part is not included, since it introduce only a 6s's factor. 

We preserve the same notations as in section 6. stands for the integral calculated from a to /3 limits; \nz) is 

the Zi,^ (z) wave function; \npmnz) is the total '^ra^n'} wave function of the basis. 

The following notations will be available from now on: 



/(n' , np, m', m, m! nig, nip, t) 



T{n'^ + l)r(n, + 1) 
r(n' + m')T{np + nig) 



/CSO ^ 
e-a;'"'>4Tl(a;)Ll7l(x)dx (82) 



which is computed by Gauss-Laguerre nummerical integration procedure, and: 



G{n'np,m',m,t) 



r(< 



T{n'p + l)Tinp + l) 
^\ni'\ + l)r{np + \m 



1) 



2ait' 



2e 



-Im'l 



(83) 



where: 



where loi is the corresponding frequency for the sphere with the same volume. Relations between ellipsoids frequencies 
and the one for the sphere can easily be obtained from volume conservation ab"^ = i?f , when a, h are the semiaxes of 
the ellipse and Ri the radius of the sphere having the same volume. So the deformation shape dependence is contained 
in ai factors. Also: 

WiP.m'a (85) 

means that the integral is performed between the limits a and h. 



(i|/(p,^)|j)(z = c) (86) 

means the integrand is taken at the variable value z = c. This last kind of expressions intervene when the spin-orbit 
differentiation operators act on the Heavyside 9{z) function, resulting the 5{z) function. 
We have for AFi : 



(n^mVjAF/^Vpmn^) = 

~[(2np + |m| + l)5n'^n, - yjnp{np + \ni\)5n^,n,-i - \J {rip + l){np + |m| + l)5n^,n,+i] 
+ \[imoujl^ - mowf)«|(z + zif\n^Yj:^ + moa;^^«|(z + zif\n^)f^^+^^ 
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—I{n'p, Up, \m\, |m|, \m\ + 1, |m| + 1, \m\ + 1, a-i^p^{z)) 



a 



-TOoW2("zl(^ - ^2) |nz)o 

For the minus terms in AVi: 



I{n'p, Up, |m|, |m|, |m| + 1, |m| + 1, |m|, aip^(2)) 



(87) 



{n'pm'n'^\AVi{zi pf)\npmn^) = 

1 
2 



^mow^^ {n'^\nz)li\[I{n'p, Up, |m|, |m|, |m| + 1, \m\ + 1, |m| + 1, alpf{z)) 



at 



-I{n'p,np, \m\, \m\, \m\ + 1, \m\ + 1, \m\ + l,afp^(^;))] 

+ ^moujl^{n'J{z + zif\n:,)z([I{n'p,np, |m|, |m|, |m| + 1, |m| + 1, |m|, afp?(2;)) 
-I{n'p,np, |m|, |m|, |m| + 1, |m| + 1, |to|, aiP^(2;))] 
Matrix elements for AV2 reads: 



(88) 



{n'pm'n'^\/^V2^\npmnz) = 

l(mo< - moa;?)(«|n,)°,^ + «|n.)r+''"*) 



~[(2np + |m| + l)5n'^n, - \Jnp{np + \m\)8n^,n,-i - \J [rip + l){np + |m| + l)J„^„^+i 



-/(n^, np, |m|, |m|, |m| + 1, |m| + 1, |m| + 1, aip„(2))] 



■mow2(n;|(z + zi)V,)° 



-{niouj^^ - moW2)«|(2 - 22) |n^)o' 
•[(^n^p - I{n'p,np, |m|, |m|, |m| + 1, \m\ + 1, \m\,alpl^{z))] 

and again for the minus terms in AV2: 



(89) 



{n'pm'n'JAV2{zi ^ Zf;pi^ pf)\npmnz) = 



^moWp, «|n^)^f ^ [/(n;;, Up, |to|, |to|, |m| + 1, |m| + 1, |m| + 1, ajp^iz)) 



-I{n'p,np, \m\, \m\, \m\ + 1, \m\ + 1, \m\ + l,alp}{z))] 

+ \^oi^l^{K\{z - Z2f\n^)li [I{n'p,np, |m|, |m|, |m| + 1, |m| + 1, \m\,a\p1{z)) 
-I{n'p,np, \m\, \m\, \m\ + 1, \m\ + 1, \m\,alp'f{z))] 

And finally for the neck potential operators, the expressions for the matrix elements are: 



(90) 



{n'pm'n'^\Vgi{zi ^ Zf;pi ^ pf)\npmnz) = 

- ^pI^ {n'Jn.)l( - \m,u^lK\{z - z,f\n.)V. 
\I{ri!p, Up, \m\, \m\, \m\ + 1, |m| + 1, \m\,alp'^{z)) 

-]-moJl{n'^\n^)li | \\I{n'p,np, \m\, |m|, \m\ + 1, |m| + 1, |m| + l,a\p1{z)) 

-I{n'p,np, |m|, |m|, |m| + 1, |to| + 1, |m| + l,alp}{z))] 



- — [I{n Up, |m|, |m|, |m| + 1, |m| + 1, \m\ + \, a^p- (z)) 
-/(n^, np, |m|, |m|, |m| + 1, |to| + 1, |m| + aiPj(z))] 



and: 



{n'pm'n'^\Vg2\npmn^) = 

yo{{n'z\nz)ll^[I{n'p,np, \m\, \m\, \m\ + 1, |m| + 1, \m\,alpl(z)) 
-I{n'p, Up, |m|, |m|, |m| + 1, |m| + 1, |m|, aiPg(z))] 
+ ("''J"-^)zf [■fl'T-p-^-p. |m|, |m| + 1, |to| + 1, \m\,alpl(z)) 
~I{n'p, Up, |m|, |m|, |m| + 1, \m\ + 1, |m|, Q;iPg(z))]} 

For the z-dependent terms of the matrix elements, the foUowing abreviations are available 

/•OO 











'Co 

where the general formula for computing these kind of intergals is: 



1^1 — 1^2 



and 7ii/(C) is the Hermite function. 



— ^y'2^'^2'^\jiy'2,U2,a2(z2-Zi) ~ il/!^,U2,a2(z2 — Z;)\ 1 ^i: Zf ^ 



{n'^\{z + ZiY\nz) zi — C^'^C^^-^ ^ -^iu[+lMi + l,ai{zi+Zf) ~ iv[+l,ui+l,ai(zi+z^)\ 



If' ■ 1 

~^~^\3v^+l,vx — \,ax{zx+Zj) ~ JiyJ+l, 1^1-1, ai(2i+Zi)J 

~r 2 Wi'I— 1,1^1 + 1. ctl(zi+z/) Jv[-l,:yi + l,ai{zi+Zi)\ 
+ ^li^lb'i/{-l, 1^1-1, Qi(zi+z/) ~ ju[-l,ui-l,ai(zi+Zi)\ 



{n'^\{z — Z2Y\nz)zi — Cv'^C^,-^ ^ -^[j i^'^ + l, 1^2 + 1 ■a2(z2- z^) ~ j}j'^ + l,V2 + l,a2(z2~Zf)\ 

'~^\iv'^ + l,V2~l,a2{z2~'Zi) ju^ + lM2 — l,a2{z2—Zf)] 



v. 



+ 1^2 ^^2 [j'l/^ —1,1/2 —1,02(22— Zi) Jfj — 1,1/2 — 1,02(22 — Z/ 
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(99) 



+ i^2[jl/^, 1/2 — 1,02(22— Zi) ~ j'^'2,'^1—'^,OC2(z2—Zf)\ } 

For [z — z^) dependent terms, one can use: 



(100) 



«|(z - Zif\n,)li = {n',\{z - Z2f\n,Yzi + 2{z2 - z^){n',\z - z^^^ + (z^ - z^f {v!,\n,)l{ 



(101) 



d_ 
'd'z 



riz ) 

I Zi 



2{jv[,v\+\,a\(z\+z f) ~ jVj, 1/1+1, ai(zi+Zi)) 
+ ^^l(ii/(, 1/1-1, Qi(2i+z/) ~ ii/;, 1^1-1, ai(zi+2i)) ^ii^/ < 

~ 2 (•^i'2'''2 + l,a;2(z2— Zi) ~ Ji/2,i/'2+l,a2(z2— z/)) 
+«^2(ii/^,i/2-l,a2(z2-Zi) -ii/^,i/2-l,a2(z2-z/)) ^i'^/ ^ 



(102) 



SPIN-ORBIT AND Q.^ MATRIX ELEMENTS 

The necessary expressions to calculate the matrix elements for spin-orbit and fi^ terms are given in this appendix. 
For the (t;i) volume, we use the abreviation: 



«;i(t;i) 



ai 





d 






d 




dz 


riz) + 
1 —00 




'd'z 



zom+a^ 

^'\^z\l^z){z = Zom + am) - {n'J\nz){z = Zim) 



(103) 



■u;2(i^i) = {.\Jnp + \m\ + ISn'^up + y/n^5n'^np-i) ■ 2mot^^i 
•(«|z + ziln,)!.™ + «|z + zi\n,)%^+,J 



ai 



+{\Jnp+ \m\ + l5n'^n, - ^/rrp5„'^n,-i) ■ 



nioUJ 



pi 





d 






d 














dz 


/ —00 




dz 



ZOm+ar, 



+ ("zl)"'z(^ = ZOm + am) - {K\nz)iz = Zim)] 



(104) 



Then we have: 
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{n'pm'n'J{n+{vi), [1 - 9{z - Zim)} + {n+{vi), [z - {zom + am)]}|npmn^) 
= Sm'm+l -wiivi) , m < 

= 5m'm+l ■ {-W2{vi)) ,171 >0 



(105) 



{n'pm'n'J{n-{vi), [1 - 6{z - Zim)} + {^+{vi), [z - (^om + a.m)]}\npmn^) 

= 5m'm^l-W2{vi) ,m<0 

= 6m'm-l ■ i-wi{vi)) ,m>0 



Because of the symmetry we have: 



(106) 



('|{fi+,/(p,z)}|)(m',m) = {\{n-,f{p,zm{m,m') 



(107) 



So, only one half of the fJs and matrix elements must be calculated. Further on, only the ri"*" matrix elements 
are given, the Cl~ part resulting from symmetry. 

{n'pm'n'J{Cl+{vi), [e{z - Zim) - 0[z - {zom + am)]]0{p - pm{z))}\npmn:,) = 
Sm'm+i{-I{n',np, \m\ - 1, |m|, |m|, |m| + 1, \m\,alp'^{z)) 



4. 



ai 



(2(< 



ZOm+ar, 



+ {n'Jnz){z = Zim) - {n'^\nz){z = zom + am) 



■[2{np + \m\)I{n'p,np, \m\ - 1, \m\ - 1, |m|, |m| + 1, |m| - l,alpl^{z)) 
+-^^Giri'p,np,\m\ - l,\m\,a\pl,{z))]} ,m < 



(108) 
(109) 



and: 



(n^m'n^|{n+(ui), [0{z - Zim) - 0[z - {zotu + am)]]d{p - Pm{z))}\npmnz) = 
Sm'm+i{I{n'p, Up, \m\ + 1, |m|, |m| + 2, \m\ + 1, |m| + 1, a^p^iz)) 



mou) 



pi 



ai 



(2(< 



d_ 

dz 



nzj^ + (n^|n^)(2 = Ztm) - {n^\n;,){z = Zom + am)j 

ZOm+Clm 



-2moUjl^ai{n'Jz + zi|n^) 
■[21{n'p, Up, \m\ + 1, |m| + 1, \m\ + 2, \m\ + 1, |m| + 1, a?p^(z)) 

^G(n^, Up, \m\ + 1, \m\,alpl^{z))]} , m > 



(110) 



(n>'<|{fi+(.;i), Mz - Zi) - e{z - zf)mp - pi{z)) - e{p - pf{z))]}\npmnz) = 
Sm'm+iil-Iin'p.np, \m\ - 1, |m|, |m|, |m| + 1, \m\,a\p1{z) 
+I{n'p, Up, \m\ - 1, |m|, |m|, |m| + 1, |m|, a?p^(2))] 

• 2mQ(Jl^ai{n'^\z + zi\)l{ 



+ 



2(< 



+ {'n''z\nz){z = Zi) - {n'^\nz){z = Zf) 
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■[2{np + \m\)[I{n'p,np, \m\ - 1, |m| - 1, |m|, \m\ + 1, |m| - l,a\pl{z)) 
-I{n'p,np, \m\ - 1, \m\ - 1, |m|, |m| + 1, |m| - 

+^[G(n^,np, |m| - l,\m\,alpl{z)) - G{n'p,np,\m\ - l,\m\,alpj{z))]]} 

,m<0 (111) 

and the same term for m > : 



{n'pm'n',\{n+ivi), [e{z - Zi) - 6{z - Zf)][e{p - pi{z)) - 6{p - p}{z))]}\npmn,) = 
Sm'm+i{[I{n'p, Up, \m\ + 1, |m|, \m\ + 2, |m| + 1, |m| + 1, a\p}{z)) 
-I{n'p,np, \m\ + 1, |m|, |m| + 2, |m| + 1, |m| + l,alp'f{z))] 

2mou)'^^ai{n'Jz + zi\n^)li 



tti 



-2TOoW^jai(n'^|z + 2i|n2)2f 

•[2(J(n;„ Up, |m| + 1, |m| + 1, |m| + 2, |m| + 1, |m| + 1, alp1{z)) 
-I{n'p, Up, \m\ + I, \m\ + 1, \m\ + 2, |m| + 1, |m| + 1, alp){z)) 

-^(G(np, rip, |m| + 1, \m\,alpi{z)) - G(n'p, Up, \m\ + 1, \m\,alpj{z))]} 

,m > 



(112) 



{n'pm'n'jn^{vi)\npmn;,) = moWp^m^^ 



(113) 



We now present the (^2) -dependent necessary terms for ns(w2) matrix elements construction: 



{n'pm'n'J{n+{v2), [0{z - Zim) - 6[z - {zQm + a.m)W - ^(p - Pm{z))]}\npmn;,) = 
Sm'm+i[^rnou)l^ai{n'Jz - Z2\n:,)l°2'^''"' {^Jrip + |m|(5„/^„p + i/np + ISn'^n^+i) 

a^'^^V^. +{n'z\nz){z = Zim)-{n'^\n^){z = ZQm+am)) 

■{\Jnp + |m|(5„^„^ - ^Jnp + 15n'^n,+l) 

+I{n'p, Up, \m\ - 1, \m\, \m\, \m\ + 1, \m\,alp'^{z)) 



ai 



d_ 

dz 



riz } + {n'Jnz){z = Zim) - {n'Jnz}{z = zo 



)) 



-2moUJ^,^ai{n'Jz - Z2\n^)l"^^+''- 

■[2{np + \m\)I{n'p, Up, \m\ - 1, \m\ - 1, |m|, |m| + 1, |m| - 1, ajpl^iz)) 
+-^G{n'p,np,\m\ - l,\m\,alp'^{z))]j ,m < 



(114) 



{n'pm'n'^\{n+{v2), [0{z - Zim) - 0[z - (^om + am)]][l - 0{p - pm{z))]}\npmn^) = 

Sm'm+i{ - 2mo(jjl^ai{n'^\z - Z2\n^)l°r^°'"' {^J Up + \m\ + l(5„^„p + V^'^n^p-i) 
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ai 



d_ 

dz 



+ {n'z\nz){z = Zim) - {n'Jnz){z = zom + a™)) 

■{^Jup + \m\ + ISn'^n, - y^Sn'^n,-l) 

-I{n'p, Up, \m\ + 1, |m|, |m| + 2, |m| + 1, |m| + 1, alp^iz)) 
2moujl^ai{n',\z - Z2\nzYzlZ^''-^ 



mouj 



P2 



ai 



•[2/(np, Up, \m\ + 1, |m| + 1, |m| + 2, |m| + 1, |m| + 1, Q!ip^(2;)) 
-^G(np, rip, |m| + 1, |TO|,aip^(z)) | ,m>0 



{n'pm'n',\{n+{v2), [e{z - Zi) - 0{z - Zf)][9{p - Mz))e{p - pf{z))]}\npmnz) = 
Sm'm+i^[-Iin'p,np, \m\ - 1, \m\, \m\, \m\ + 1, \m\,alpi{z)) 
+I{n'p, Up, \m\ - 1, \m\, \m\, \m\ + 1, \m\,alp}{z))] 
■ 2moujl^ai{n'J\z - Z2\nz)l( 



P2 



ai 



d_ 

dz 



+ {n'z\nz){z = Zi) - {n'^\nz){z = Zf)^ 
+2moU!l^ai{n'Jz - Z2\nz)zi 

■[2{np + \m\){I{n'p, Up, \m\ - 1, |m| - 1, |m|, |m| + 1, |m| - 1, alp^{z)) 



-I{n'p, Up, |m| - 1, |m| - 1, |m|, |m| + 1, |m| - 1, alpj{z))) 



+ 



2ai 



{G{n'p,np, \m\ - 1, \m\,alp'^{z)) - G{np,np, \m\ - 1, \m\,alp'f{z)))]j 



,m < 



(115) 



(116) 



{n'pm'n'Mf^+{v2), [e{z - z,) - e{z - Zf)][9ip - p^iz))9ip - pf{z))]}\npmnz) 
Sm'm+i^[I{'n''p,np, \m\ + 1, \m\, \m\ + 2, \m\ + 1, \m\ + \,a\p1{z)) 
-I{n'p, Hp, \m\ + 1, |m|, |m| + 2, |m| + 2, |m| + 1, |m| + 1, alp}{z))] 
■ 2mou)l^ai{n'^\z - Z2\nz)li 



ai 



-(2(< 



d_ 

dz 



' + {n'z\nz){z = Zi) - {n'Jnz){z = zj)^ 
-2moUjl^ai{n'^\z - Z2\nz)li 

■[2{I{n'p, Up, \m\ + 1, |m| + 1, |m| + 2, |m| + 1, |m| + 1, alp1{z)) 
-I{n'p, Up, |m| + 1, |m| + 1, |m| + 2, |rn| + 1, |m| + 1, alp^f{z))) 

-^(G(n^,np, + 1, \m\,alpf{z)) - G{n'p,np, \m\ + 1, \m\,alpj{z)))]^ 

,m > 



The diagonal term generated by flz is: 



(117) 



{n'pm'n'Jftz{v2)\npmnz) = moujl^rndm'm 



(118) 



For the neck region elements we have: 
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(nX<|{n+Ki), [e{z - Zi) - 6{z - Zf)][e{p - pi{z)) - 9{p - p/(^))]}|n,mn,) = 
5m'm+imouig{ { — [I{n'np, \m\ - 1, |m|, |m|, |m| + 1, \m\,alpi{z)) 



-I{n'p,np, \m\ - 1, |m|, |m|, |m| + 1, \m\,aj^Pf{z))] 
-P3[l{n'p, Hp, \m\ - 1, |m|, |m|, |m| + 1, |m| - ^, aiP?(2)) 

-7(n^,np,|m| - l,|m|, |m|, |m| + 1, |m| - ^,alpj{z))]} 



2 n' 



d_ 

dz 



{n'^\n,){z = Zi) - {n'^\n^){z = Zf) 

-2{ai{2(np + \m\)[I{n'p, Up, \m\ - 1, |m| - 1, \m\, \m\ + 1, |m| - 1, a?p,2(^)) 
-I{n'p, Up, \m\ - 1, \m\ - 1, |m|, |m| + 1, |m| - 1, alp}{z))] 
-[I{n'p,np, \m\ - 1, |m|, |m|, |m| + 1, \m\,alp'^{z)) 
-I{n'p, Up, \m\ - 1, |m|, |m|, |m| + 1, \m\,alp}{z))]} 

+ 7:[G{n'p,np, \m\ - 1, \m\,alpf{z)) - G{np,np, \m\ - 1, \mlalp){z))] 



■{n'^\z - Z2,\nz)li}^ ,m 



< 



(119) 



{n'yn',\{n+{vgi), [e{z - Zi) - e{z - Zf)][e{p - pi{z)) - e{p - pf{z))]}\npmn,) 
5m.'m.+i'mQUj1^[I{n'p, Up, \m\ + 1, |m|, |m| + 2, |m| + 1, \m\ + 1, a\p1{z)) 
-I{n'p,np, \m\ + 1, |m|, |m| + 2, |m| + 1, |m| + l,alp^f{z))] 





9 








9z 



-2a(<| 
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d 








dz 



-P3,[l{n'p, Up, |m| + 1, |m|, |m| + 2, |m| + 1, \m\ + -, aiPj (2;)) 
J(n|,, Hp, \m\ + 1, |m|, |m| + 2, |m| + 1, |m| + ^, aip^(z))] 

Zf 

+ Kh^X^ = - {n'z\nz){z = Zf) 
+2{2ai[/(n^, rip, |m| + 1, |m| + 1, |m| + 2, |m| + 1, |m| + 1, Q;?p,2(0)) 
-/(n;,, rip, \m\ + 1, |m| + 1, |m| + 2, |to| + 1, |m| + 1, a\p}{z))] 

-\[G{n'p, Up, \m\ + 1, \m\,alpl{z)) - G{n'p, Up, |m| + 1, \m\,alp}{z))]} 
■{n'Jz - Z3\nz)t(^ ,m>0 



(120) 



{n'pm'n'^\{n;:{vgi), [e{z - Zi) - ^(2; - 2;/)][6'(p - pi{z)) - e{p - pf{z))]}\npmn;,) = 
moujlmS„i'jn{[I{n'p, Up, \m\, \m\, \m\ + 1, |m| + 1, \m\,alp'^{z)) 
-I(n'p, Up, |to|, |to|, \m\ + 1, |to| + 1, \m\,alp}{z))] 

-P3ai[l{n'p,np, \m\, \m\, \m\ + 1, |m| + 1, \m\ - ^,alpi{z)) 

-I{n'p,np, \m\, \m\, \m\ + 1, \m\ + 1, \m\ - ^,alpj{z))]} (121) 



23 



* rgherg@ifin.nipne.rq 
[1] E. Merzbacher, Quantum Mechanics (John Wiley and Sons Inc., New York, 1961). 
[2] J.Mamhn and W. Greiner, Z. Phys. 251, 431 (1972). 
[3] W. vonOertzen, Z. Phys. 354 1, 37 (1996). 
[4] A. Thiel, J. Phys. G: Nucl. Part. Phys. 16, 867 (1990). 
[5] J. W. Park, W. Scheid and W. Greiner, Phys. Rev. C21, 958 (1980). 

[6] D. N. Poenaru, M. Ivascu and W. Greiner, Fission and Beta-Delayed Decay Modes, vol. Ill, (CRC Press, Boca Raton, 
Florida, (1989); D. N. Poenaru and W. Greiner, Nuclear Decay Modes (Institute of Physics Publishing, Bristol, England, 
ch.6, (1996)) 

[7] R. A. Gherghescu, J. Skalski, Z. Patyk and A. Sobiczewski, Nucl. Phys. A651, 237 (1999). 

[8] W. I. van Rij and C. T. Hess, Nucl. Phys. A142, 72 (1970). 

[9] E. Badralexe, M. Rizea and A. Sandulescu, Rev. Roum. Phys. 19, 63 (1974). 
[10] P. Moeller, J. R. Nix and K.-L. Kratz, Atomic Data Nucl. Data Tables 66, 131, (1997). 
[11] A. C. Wahl, Atomic Data Nucl. Data Tables, 39, 1, (1998). 



24 



FIGURE CAPTIONS 



Figure 1 Nuclear shape and deformation parameters for necked-in intersected ellipsoids. Oi, O2 and O3 are the 
centers of the two deformed fragments and of the neck sphere respectively. The four geometrical parameters which 
vary arc the two ratios of the ellipsoid semiaxes, bi/ai and 62/^2, the neck sphere radius Rs and the distance between 
the fragment centers R. 

Figure 2 Sequences of shapes when geometrical parameters are varied. Neck sphere radius vary from small (upper 
row - close to fusion shapes) through intermediary (second and third row) up to large values for very elongated shapes. 
Along one row, the distance between centers is increased for the same value of the neck. Semiaxes ratios and mass 
asymmetry here are kept constant. 

Figure 3 Matching potential ellipsoids: MPE between the two fragments (shaded area) and between the neck potential 
and the fragment potentials (grey area). The influence of each of the potentials is emphasized in the upper part. The 
same figure with the surface functions for the regions of interest (fragments: pi{z) and P2{z)', neck: Pg{z); MPE: p{z); 
MPEl: pmij{z), Pmis{z), MPE2: pm2s{z) and Pm2j{z) and the 2;-values used in calculations are marked. 

Figure 4 The evolution the two fragment potentials Vi{p,z), V2{p,z) and the neck potential Vg{p,z) for the p- value 
at the two tangent points of the neck with the fragment ellipsoids: p{zci) - upper figure and p{zc2) when the distance 
between centers increases. One can observe how the three potentials are tangent all along the splitting. 

Figure 5 Variation of the heavy fragment potential Vi{p,z) and neck potential Vg{p,z) along the O1O3 direction. 
One can observe that inside the neck sphere - hence inside MPEl and MPE2, the neck potential is higher, acting as 
a pressure element on the nuclear shape. The two crossing points at 2i?3 distance on O1O3 direction mark the limits 
of MPEl along the O1O3 direction. 

Figure 6 Shell corrections calculated for the levels presented in Fig. 6, as a function of the distance between centers 
R. The ground state is practically common. A first bump appears in all cases; then for larger R3 there is a decrease 
around i?=10 fm. The second bump, which usualy generates a second maximum in the fission barrier, lies around 15 
fm, but is lower for small R3. 

Figure 7 Two level spectra with different mass asymmetry from the same superheavy parent ■'''^122 are presented on 
the upper part, for i?3=4 fm. The lower part of the figure shows the corresponding shell correction energy EsUeii as a 
function of R. First bump is present in both cases, eventhough a little shifted. At the end of the splitting, asymmetric 
channel seems more favorable from microscopic point of view, as shell correction energy is some 4 MeV lower. 

Figure 8 The level schemes for two fission channels from ^'''^Cf are displayed in the upper part of the figure. Levels 
are almost likely in the first part of the splitting as can be deduced from the corresponding shell corrections. When the 
two fragments continue to split, the individual wells influence results in a rather large difference in shell corrections, 
beyond R=7 fm. 
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Figure 4 
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